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Abstract. Two-pion Hanbury-Brown-Twiss (HBT) correlations for p+p and central Pb+Pb 
collisions at the Large-Hadron-Collider (LHC) energies are investigated with the ultra- 
pH | relativistic quantum molecular dynamics model combined with a correlation afterburner. The 

transverse momentum dependence of the Pratt-Bertsch HBT radii Riong, Rout, and R S ide is 
extracted from a three-dimensional Gaussian fit to the correlator in the longitudinal co-moving 
system. In the p+p case, the dependence of correlations on the charged particle multiplicity 
and formation time is explored and the data allows to constrain the formation time in the string 
fragmentation to 77 < 0.8 fm/c. In the Pb+Pb case, it is found that R ou t is overpredicted 
by nearly 50%. The LHC results are also compared to data from the STAR experiment at 
RHIC. For both energies we find that the calculated Rout/ Raids ratio is always larger than 
data, indicating that the emission in the model is less explosive than observed in the data. 



1. Introduction 

The properties of strongly interacting matter are described by the theory of Quantum-Chromo- 
Dynamics (QCD). To explore the details of QCD matter under extreme conditions, one needs 
to compress and heat up QCD matter to regimes present microseconds after the Big Bang. 
Today these conditions can only be found in the interior of neutron stars or created in heavy-ion 
collisions at relativistic energies. Over the last decade the experimental programs at the Super 
Proton Synchrotron (SPS) and at the Relativistic Heavy Ion Collider (RHIC) have provided 
exciting pioneering data on the equation of state (EoS), the transport properties of the matter 
created and its spatial distributions [Xl [21 El SI El El [3 [HJ [91 [10]. In addition, the Large 
Hadron Collider (LHC) at CERN had been designed, installed, tested, and repaired in the 
past two decades and finally, started normal operation in the end of the year 2009. Since 
then, a tremendous amount of experimental data in various aspects of high energy physics 
has been obtained and received much attention by theoretical physicists. And, the extracted 
bulk properties of the high temperature fireball created in such ultra-relativistic collisions have 
provided unprecedented information for fundamental investigations of the phase diagram of 
QCD. Here we want to explore the expansion properties of the created matter by investigating 
the spatial shape of the fireball. Although it is known that one can not measure the emission 



time pattern and/or the spatial profile of the source directly, a well-established technique, called 
"femtoscopy" or "HBT" (see e.g. and references therein) can be employed to obtain this 
information. Femtoscopy has been extensively used in the heavy ion community since it provides 
the most direct link to the lifetime and size of nuclear reactions. The ALICE collaboration has 
published first results of two-pion Bose-Einstein correlations in both p-p [12] and central Pb-Pb 
|13| collisions at LHC energies in the beginning of the year 2011. These experimental results 
have attracted the research interest of several theoretical groups [1'6\ 114 4 IT~5 ^, 116]. whose models 
are based on hydrodynamic/hydrokinetic and hadronic microscopic approaches. 

It is often argued that the particle emitting system in p+p collisions is too small to create 
a medium that exhibits bulk properties, this is different at a center of mass energy of \fs= 7 
TeV. Here, the particle multiplicity is about the same as in nucleus-nucleus collisions, studied 
at RHIC. These data suggest that space momentum correlations are developed even in p+p 
collisions as soon as high particle multiplicities are achieved. Thus it is worthwhile studying 
the dependence of HBT observables on the event multiplicity. As the system created in p+p 
collisions is still small, an essential quantity that influences the particle freezeout radii is the 
formation time in flux tube fragmentation. The formation time sets the scale for a minimum 
value of the source lifetime - of course followed by resonance decay and rescattering. In this 
paper we suggest that the recent LHC data on p+p collisions allow us to determine the formation 
time in the flux tube break-up. 

In this paper we show results for the HBT radii of two-pion correlations from p+p and 
central (< 5% of the total cross section or) Pb+Pb collisions at LHC energies ^/snn = 7 and 
2.76 TeV, respectively, from the ultra-relativistic quantum molecular dynamics (UrQMD) model 
|17[ \18\ \19\ [20] . The calculations are compared to ALICE data as well as to those at the RHIC 
energy y/sJm = 200 GeV. 

The paper is arranged as follows: In Section 2, a brief description of the UrQMD model and 
the treatment of the HBT correlations as well as the corresponding Gaussian fitting procedure 
are shown. Section 3 gives the main results of the model calculations for p+p and Pb+Pb 
collisions at LHC energies. Finally, in Section 4, a summary is given. 

2. Brief description of the UrQMD model and the HBT Gaussian fitting procedure 

UrQMD \19\ [20"1 [2T] is a microscopic many-body approach to p+p, p+A, and A+A interactions 
at energies ranging from SIS up to LHC. It is based on the covariant propagation of mesons and 
baryons. Furthermore, it includes rescattering of particles, the excitation and fragmentation of 
color strings, and the formation and decay of hadronic resonances. At LHC, the inclusion of 
hard partonic interactions in the initial stage is important and is treated via the PYTHIA [22] 
model. 

Formation time denotes the time it takes for a hadron to be produced from a fragmenting 
string. A very common model to describe such a flux tube fragmentation is the Lund string 
model [23]. In the Lund model the formation time consists of the time it takes to produce an 
quark-antiquark pair and the time it takes for that pair to form a hadron. For the Lund model, 
both of these times are proportional to the transverse mass of the created hadron and inversely 
proportional to the string tension. For simplicity, UrQMD uses a constant formation time of 
Tf = 0.8 fm/c for hard collisions. Since HBT probes the freezeout distribution, the extracted 
radii are sensitive to the value of Tf in small systems at high collision energies where the majority 
of particles are produced from flux tube fragmentaion. 

In the present study, the cascade mode of the latest version (v3.3) of UrQMD is used (for 
details of version 3.3. see |19[ [20]). Some predictions and comparison works with data from 
reactions at LHC have already been pursued based on this version and showed encouraging 
results for the bulk properties [24^ 125] . After the UrQMD simulation, all particles with their 
phase space distributions at their respective freeze-out time (last collisions) are put into an 



analyzing program using the formalism of the well known "correlation after-burner" (CRAB) 
|26| . At last, the constructed two-pion HBT correlator (regardless of charge) in the longitudinally 
co-moving system (LCMS) [27} [28] without influence of residual interactions is fitted (using the 
\ 2 method) with a three-dimensional Gaussian form expressed as 



C(q, k) = 1 + A(k)exp 



E 



i,j=out,side,long 



(1) 



In Eq. ([T]), A represents the fraction of correlated pairs [IT] and qi is the pair relative momentum 
q = pi — p2 in the i direction, pi being the momenta of the particles, long is the longitudinal 
direction along the beam axis, out is the outward direction along the transverse component of 
the average momentum k of two particles (ky = |pit+P2t|/2) and side is the sideward direction 
perpendicular to the afore mentioned directions. The effect of cross terms with i / j on the 
HBT radii is found to be negligible in the present fits when the pseudorapidity cuts \n\ < 1.2 
and \t)\ < 0.8 are used for p+p and Pb+Pb cases, respectively, as in experiments, and is not 
discussed in this paper. 

For central collisions, the HBT radii are, except for an implicit k^ dependence, related to 
regions of homogeneity by [29] 
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Rl de = (y 2 ), (3) 

and 

Ring = ((Z ~ htf) = (z 2 ) ~ 2 (fotz) + (f5 2 L t 2 ) , (4) 

where x, y, z and t are the spatio-temporal separation of the particles in a pair and (3= k/feo- 
If no space-momentum correlations are present the regions of homogeneity and the source size 
coincide. In central collisions the relation (x 2 ) ~ (y 2 ) is satisfied. Thus R 2 ut and R 2 ide differ 
mainly in the last two terms of Eq. ([2|). The first of these two terms is dependent on the 
strength of the correlation of emission time and transverse emission position, while the second 
one is especially sensitive to the particle emission duration. 

For the p+p case, the analysis is done differentially for the different event multiplicities listed 
in table [1] Generally, the average dN c h/dn from UrQMD is 15% smaller than the one measured 
by the ALICE collaboration, in the same charged multiplicity classes. 



3. Results 

3.1. p+p collisions at yJspjN = 7 TeV 

In this subsection, the results on HBT radii in p+p collisions from the UrQMD model are 
calculated and compared to ALICE data [12]. In Fig. [1] the projections in the out, side 
and long directions of the 3D correlation function together with a projection of the fit for 
kx = |ky| = 0.3 — 0.4 GeV and N^^ g 2 d = 23-29 are shown as an example. It can be seen 
that the calculated correlation function (shown as dots) is well described by a Gaussian fit 
(lines). However, oscillations of the correlation function are present at larger q. This indicates 
that there is a non-gaussian component in the underlying separation distribution of the pion 
freezeout points. Here we will not investigate this interesting question further. 

The kx dependence of the HBT radii extracted from the UrQMD freezeout distribution is 
presented in Fig. [5] for the dN c h / dn classes defined in table [T] in the pseudorapidity interval 
\n\ < 1.2 in comparison to the ALICE data. The UrQMD calculations are presented for dif- 
ferent values of the formation time tj (t/= 0.3 fm/c, dashed lines; ry= 0.8 fm/c, full lines; 



Table 1. Table of the investigated multiplicity intervals. The first column shows the interval 
boundaries, the second column the mean charged particle multiplicity per unit of pseudorapidity 
(dN c h j drj) from events with at least one charged particle in \rj\ < 1.2 from UrQMD. The third 
column shows the same quantity from ALICE data |12j . 



N \n\<i-2 

charged 


UrQMD {dN ch /dr,) H<l 2 


ALICE {dN ch /dr,) [vl<h2 


1-11 


2.52 


3.2 


12-16 


5.74 


7.4 


17-22 


8.01 


10.4 


23-29 


10.72 


13.6 


30-36 


13.65 


17.1 


37-44 


16.74 


20.2 


45-57 


20.94 


24.2 


58-149 


27.57 


31.1 




Figure 1. The dots represent projections of the three-dimensional correlation function for 
/c T =0.3-0.4 GeV and N^^ d = 23-29. The lines represent a x 2 fit of Eq. (JD to the correlation 
function. Both the result of the fit and the correlation function are integrated over a range of 
qi = ±0.17 GeV in the other directions for the purpose of projection. 



tj=2 fm/c, dotted lines). For t/=0.3 fm/c one obtains a good description for R ou t, while R s id e 
is slightly under predicted and the values for Ri ong are in line with data from ALICE. The 
choice Ty=0.8 fm/c leads to a slight overestimation of Rout] however, it leads to a reasonable 
description of R s ide data. Also the kx behaviour in R[ ong is much closer to the behaviour of the 
data. In contrast, a formation time of T/=2 fm/c, leads to a drastic overestimation of the data 
for all radii. Although there are discrepancies between model and data for all values of rj, the 



sensitivity to Tf is much larger than those discrepancies. Therefore, the present ALICE data 
allows to constrain the formation time to values of rt ~ 0.3-0.8 fm/c. 




Figure 2. The lines represent HBT radii in p+p collisions at y / s= 7 TeV from UrQMD for 
different multiplicities and formation times. The various line styles refer to results for tj = 0.3 
fm/c (dashed), tj = 0.8 fm/c (default full lines) and tj = 2 fm/c (dotted). The colors represent 
the multiplicity classes as defined in table [TJ The points are data from the ALICE experiment 

The overall shape of the radii as a function of multiplicity and kx is also interesting. The 
Rside radii (see Fig. [21 middle) from UrQMD and in the data are flat as a function of kx for low 
multiplicity events. With increasing multiplicity, the radii develop a decrease towards higher 
kx- This is exactly the behavior one would expect for the development of space-momentum 
correlations with rising event multiplicity [16]. For R ou t (Fig. [21 left) however, there is a 
kx dependence present in all multiplicity classes. Thus, only the development of radial flow 
with rising particle multiplicity seems insufficient to explain the kt dependence. Since R ou t 
and Ri on g contain lifetime contributions of the source and R s %de does not, there seem to be an 
additional non-trivial kt and multiplicity dependence in the emission duration needed to explain 
the difference in the behavior of R ou t, Rside and Riong- This additional correlation might be due 
to an additional momentum dependence in Tf apart from the trivial Lorentz boost. This would 
lead to a direct effect on the emission duration, because it changes the particles production 
spacetime points. 

3.2. Pb+Pb collisions at ^ajfij = 2.76 TeV 

Firstly, the correlation functions are also calculated in bins of the transverse momentum kx- 
Fig. [3] shows the projections of the three-dimensional correlation function (dots) and of the 
respective fit (lines) for the kx bin 0.2 — 0.3 GeV/c. It is seen clearly that the correlator in 




Figure 3. Projections of the three-dimensional correlation function (dots) and of the respective 
fit (lines) for the hr bin 0.2 — 0.3 GeV/c and |t?| < 0.8. When projecting on one axis the other 
two components are restricted to the range (-30 < q < 30) MeV/c. 

sideward direction can be described by a Gaussian form fairly well. However, it deviates slightly 
from a Gaussian in the other two directions, especially in the longitudinal direction, as found and 
discussed in previous publications for HICs at lower energies [30]. At LHC, the fraction of excited 
unstable particles is much larger than at lower beam energies, therefore, the non-Gaussian effect 
is more pronounced in the current calculations. At RHIC energies, the non-Gaussian effect 
was also seen in the experimental HBT correlator, especially in the longitudinal direction pQ. 
Furthermore, when comparing our fitting result in Fig. [3] with that observed in experiment (in 
Fig. 1 of Ref. [13]). it is seen that the non-Gaussian effect is stronger in our calculations than in 
experiment. This might be due to the omission of all potential interactions between particles in 
the current cascade calculations. Support for this interpretation was found in Ref. [30], where 
the consideration of a mean-field potential plus Coulomb potential significantly reduced the non- 
Gaussian effect on correlators of pions from HICs at AGS energies. Therefore, both a dynamic 
treatment of the particle transport with a proper EoS for the QGP phase and the hadron phase, 
and further theoretical development of the fitting formalism are equally important for a more 
precise extraction of spatio-temporal information of the source |31| . 

Besides the non-Gaussian effect, the contribution of the correlation between the emission time 
and position to the HBT radii, especially in the outward direction, has been paid more attention 
in recent years since it closely relates to the stiffness of the EoS of nuclear matter especially at 
the early stage of the whole dynamic process [32J E31 EI] • It was found that, even in the cascade 
calculation, there exists a visibly positive correlation between the emission time and position 
[35] • And, the most important contribution to R ou t comes from the emission duration term. It 
implies that both a shorter duration time and a stronger x — t correlation lead to a smaller R ou t 
value, which will be shown in Fig. H] specialized for the result of HBT radii. 

Fig. 2] shows the kx dependence of the HBT radii R ou t, Rside-, and Ri ong , as well as the ratio 
Rout/Rside, extracted from the Gaussian fit to the two-pion correlators. The UrQMD cascade 
calculations for central Pb+Pb collisions at LHC energy -y/J/vTV = 2.76 TeV (lines with crosses) 
and central Au+Au collisions at RHIC energy s/snn = 200 GeV (lines with up-triangles) are 




Figure 4. hp dependence of pion HBT radii R out [panel (a)], R s ide [(b)], and Ri on g [(c)], as 
well as the ratio Rout/ Rside [(d)], for central (ct/ctt < 5%) Pb+Pb collisions at LHC energy 
y/SNN = 2.76 TeV. For comparison, parameters for central (a far < 5%) Au+Au collisions 
at RHIC energy y / Sivjv = 200 GeV are also shown. Lines with up-triangles and crosses are 
for model calculations while scattered symbols are for experimental data of STAR/RHIC and 
ALICE/LHC collaborations taken from [UQ2]. Lines with circles and diamonds show results 
with an artificially decreased emission duration by a factor of t sca \ e = 2 and 3, separately, in the 
analysis of correlation function. 

compared to corresponding experimental data by ALICE/LHC (open stars) and STAR/RHIC 
(open diamonds). A strong decrease of the three HBT-radii with kx is seen both in experiments 
and in the UrQMD calculations for HICs. This implies a substantial expansion of the source 
and is qualitatively captured by the UrQMD dynamics. Following experimental results, the 
calculated HBT radii for Pb+Pb at LHC are found to be larger than those for Au+Au at RHIC. 
The largest increase exists in the longitudinal direction, which is also seen by the experiments. 
Although the comparison of the calculated HBT radii Ri on g and R s ide with corresponding data 
at RHIC is fairly well, it gets worse when going to the higher LHC energy. At LHC the calculated 
R S ide values at all kr are found to be slightly smaller than data, while Riong and R ou t values 
are larger than data. Together with large calculated R ou t, the emission time related quantity 
Rout /Rside is found to be markedly larger than the data. 

From Eq. ([2]) it is clear that the value of R ou t is strongly dependent on the emission duration 
of the particles. To further investigate the contribution of the emission duration to the HBT 
radii, we artificially decrease it by rescaling the relative time t to the "effective source center 
time" t (=< U >) by t = tj — t — > t' = (tj — t)/t sca i e in the calculation of the correlation function 
at LHC energies. This effectively changes Eq. ([2]) to 

Rout = \( x ~ PTt ) > = {X } - 2- + -2 • (5) 

Iscale l scale 

The results for this calculation are presented as lines with circles (t SC aie = 2) and with diamonds 



(tscale = 3) in Fig. SI The artificially decreased emission duration leads to smaller R ou t values in 
all kx bins but leaves R s ide unchanged, as expected. Overall it results in an improved agreement 
with the data of R ou t/R s ide ratio. From Fig. H] it is also found that Riong is overestimated at 
LHC. Since Riong is mainly related to the lifetime of the source, it implies that this lifetime is 
also overestimated by UrQMD. Other calculations in |13|. [36] show that UrQMD overestimates 
the source lifetime by a factor of ~ 2 — 3 when compared to LHC data. The overestimation 
of both Rout and Ri ong can be attributed to the known fact that the pressure in the early 
stage is not strong enough in the cascade model calculations. A higher pressure would lead 
to a more explosive expansion, a stronger phase-space correlation, and a faster decoupling of 
the system, thus leading to smaller regions of homogeneity. For more discussion we refer the 
reader to [31, 33J. With the improved integrated Boltzmann + hydrodynamics hybrid approach 
(2U1 ESJ [HTJ, [35] , where various EoS of nuclear matter during the hydrodynamic evolution may be 
treated consistently and a decoupling supplemented by realistic three-dimensional hypersurfaces 
we hope to get a satisfactory solution in the near future. 

4. Summary 

To summarize, two-pion HBT correlations for p+p and central Pb+Pb collisions at the LHC 
energies are investigated with the UrQMD model combined with a correlation afterburner 
CRAB. The transverse momentum dependence of the Pratt-Bertsch HBT radii Riong, Rout, 
and R s ide is extracted from a three-dimensional Gaussian fit to the correlator in the LCMS 
system. It is seen that the correlator in sideward direction can be described by a Gaussian form 
fairly well while non-Gaussian effect is seen in both longitudinal and outward directions. 

In the p+p collisions, we have shown that the data provides rather direct access to the particle 
formation times in the flux tube fragmentation. The sensitivity to tj is large enough compared 
to the model uncertainties to find that a value ry = 0.3 — 0.8 fm/c is strongly favored compared 
to larger values for rt. Values of Tf > 2 fm/c can be ruled out from the present analysis. 

In the Pb+Pb case, both the transverse momentum k? dependence and the beam energy 
(from RHIC to LHC) dependence of the HBT radii R ou t, Rside, and Riong exhibit qualitatively 
the same behavior as found in the experiments. However, the calculated R ou t/ 'Rside ratios at all 
kx bins are found to be larger than in the data, both at RHIC & LHC. We traced this finding 
back to the explosive dynamics of the fireball at LHC which results in both a shorter emission 
duration and a stronger time-space correlation than modeled here. 
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